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Abstract 

An analytical model for the soliton-potential interaction is presented, by constructing a collective 
coordinate for the system. Most of the characters of the interaction are derived analytically while 
they are calculated by other models numerically. We find that the behaviour of the soliton is like 
a point particle 'living' under the influence of a complicated potential, that is a function of soliton 
velocity and the potential parameters. The analytic model does not have a clear prediction for the 
islands of initial velocities in which the soliton may reflect back or escape over the potential well. 

1 Introduction 

Topological solitons are widely use as models for the description of particles generated as non- 
trivial solutions of nonlinear field theories. Skyrmions are solitons which are used as a model of 
hadrons. Some solitonic solutions appear in two-dimensional Quantum chromodynamics (QCD2). 
In bosonized QCD2 these type of solutions emerge as describing baryons and quark solitons, re- 
spectively. The generalized sine-Gordon model arises as the low-energy effective action of bosonized 
QCD2 for unequal quark mass parameters. Also in the strong-coupling limit the static classical 
soliton which describes a baryon in QCD2 turns out to be ordinary sine-Gordon kink. Modeling 
of optical self focusing phenomena, magnetic fluxes in real Josephson junctions are examples from 
other branches of physics. 

In the meantime, dynamical evolution of a soliton during the interaction with potentials is an 
important phenomenon from the mathematical point of view and also because of its applications. 
Most of the researches are in base of numerical studies because such systems are generally non- 
integrable. So it is clear that we need suitable models with analytic solutions to test the validity 
of such phenomenon and predict their behaviour. 

In this paper an analytic model for the interaction of sine-Gordon solitons with defects is 
presented and the results are compared with numerical simulation outcomes from other models. So 
we need a brief review of the available models which presents in section 2. The analytic model is 
introduced and will be solved in section 3. Presented model will be compared with other models in 
section 4. The results for the soliton-barrier system are presented in sectionS. The results will be 
compared with the predictions of other models in this section too. In section 6 soliton-well system 
is discussed. Some conclusion and remarks will be presented in section 7. 

2 Models for soliton-potential systems 

model 1: The potential generally arises from medium properties. The effects of medium disorders 
and impurities can be added to the equation of motion as perturbative terms. In this method, 



scattering of a soliton by a single impurity has been modeled asfT] [2] 



dU 

(j)u-(l>xx + {l + (jSix)) — = (1) 

wheretJ denotes the strength of impurity and ^ — sincj) for the sine-Gordon model. For an 
attractive potential well,cr is negative (cr < ) and for a barrier cr is a positive number (cr > 0). 

The impurity has been added as an external potential in this model . The interaction can be 
analyzed in term of some degree of freedom for the soliton ( position of the center of the soliton) 
and an impurity mode for the external potential. 

In this approach the impurity causes the interaction of a soliton with an effective potential. In 
particular soliton can be trapped by an attractive potential because of energy loss due to radiation. 
In this model the impurity is not a rigid object. It has a localized oscillating state, so-called impurity 
mode. 

In the absence of the impurity {a — ) equation (1) has an exact one soliton solution as 

(j)k=4: arctan (^exp (^^L^£^^^ (2) 

where X(t) — Xq — Vt and V is the soliton velocity. 

If we linearize equation around its ground state, we have 

<Ptt - (jjxx + (1 + a5{x)) = (3) 

which has a localized oscillating mode 

4>impurity(.X,t) ^ a(t) CXp (^"O-Y-^ ^ ^ 

Two dynamical variables, X(t) and a(t) explain the dynamics of the soliton-potential system. 
One can describe the soliton-impurity interaction by substituting cj) — (pk + i'impurity into the 
lagrangian of the system and integrating over the variable 'x' [2]. After that, the kink coordinate 
X(t) and impurity mode a(t) are considered as collective coordinate variables and their evolution 
describe the situation of the soliton during the interaction. Therefore the soliton is changed to a 
point particle with an effective mass of me// = 8 in the effective potential V{X) — -^^^tx- This 
potential creates the effective force 

cosh-^{X) 

Model 2: The effects of the potential also can be taken into account by making some parameters 
of the equation of motion (or lagrangian) to be as functions of space or time '^Ai- In this approach, 
a finite size, finite strength potential is included by appropriately modifying the coefficient of the 
nonlinear term in the lagrangian or equation of motion. The effective lagrangian from this model 
for the sine-Gordon soliton-potential system is 

L^d^,(t)d>'(t>~X^{\-cos4>) (6) 

with solution 
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k = 4 arctan ( exp ( X^-^^^ \ \ (7) 



A is chosen as 
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where p is the width of the potentiaL For Ao < 1 we have a potential weU and Aq > 1 describes a 
potential barrier. A delta-like potential with the strength of eo is constructed with the constraint 
Xop = eo. 

Model 3: One can add such effects to the lagrangian of the system by introducing a suitable 
nontrivial metric for the back ground space-time, without missing the topological boundary con- 
ditions [SI m [7l [8] . In other words, the metric carries the information of the medium. The general 
form of the action in an arbitrary metric is: 

L{<t>,df,<t>)^/^(rxdt (9) 



where "g" is the determinant of the metric g^'^{x). Energy density of the "field + potential" can 
be found by varying " both" the field and the metric 7 . For the lagrangian of the form 

L^^-d^ct^d^4>-u{<t^) (10) 

the equation of motion becomes [3 [3] 

^ ^ ^d^4>d^ci, + 5^</)9'' V^) + = (11) 



The suitable metric in the presence of a weak potential V(x) is [SJ [HI [7] : 

»,) (12) 

The equation of motion (11) (describes by Lagrangian (10)) in the background space-time (12) is 

n ^ v( \^ ^'"^ _ ^''^ _ 1 dV{x) dc/> du{^) 

U+l/WJg^2 9.^2 2\1 + V{x)\ dx dx^ 90 " ^ ' 

For the sine-Gordon model, we have U{(t>) = 1 — coscf). A potential of the form V{x) = ae^^'-^"'^-'^ 
has been chosen in [7] while a square shape potential has been used for simulations in [8]. In the 
above potential, parameter "a" controls the strength of the potential, "b" represents its range, and 
"c" indicates the center of the potential. If a > 0, the potential shows a barrier and for a < the 
potential acts as a potential well. 



3 collective coordinate system for model 3 

The center of a soliton can be considered as a particle, if we look at this variable as a collective 
coordinate. The collective coordinate could be related to the potential by using one of the above 
models. The third model is able to give us analytic solution for the evolution of the soliton center 
during the soliton-potential interaction. 

Here we work on the sine-Gordon model with its one soliton solution of (2). By inserting the 
solution (2) in the lagrangian (10) and using metric (12), with adiabatic approximation (TJ [2] we 
have 

L ^ I {^f , , ' (14) 

2^^ ^' cosh^ {x - X{t)) cosh^ {x- X{t)) ^ ' 

For the weak potential V(x) (14) becomes 



l-h-T/(a;)^ ^^^^ (l + -V{x)\ 5— ^ (15) 

2 ^ 7 cosh^ (x - X(t)) V 2^7 cosh^ (x - X{t)) ^ ' 
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X(t) remains as a collective coordinate if we integrate (15) over variable x 



f . , -2 -2 f V(x)dx „ „ /■ V(x)dx 

J J cosh^ (x- X(t)) J 



(16) 



cosh^ (a; - X(i)) J cosh^ {x - X{t)) 

The equation of motion for the variable X(t) results from the (16) 

J cosh^{x-X{t)) V J J cosh^ {x - X{t)) ^ ' 

It is a general equation for the any kind of potential. If we take the potential V{x) = —e5{x) then 
(17) becomes 

■■ / 2,e \ (iX^ \ 4esinhX , , 

8X 1 2— + — - + 1 . = (18 

V 4cosh^X/ y 2 j cosh^X ' 

The above equation shows that the energy peak of the soliton moves under the influence of a 
complicated force which is function of its position and velocity. Note that an effective force in 
the form of equation (5) of model 1 appears in equation (18) when the soliton velocity is small 
{X ^ Q). If e > we have a barrier and e < creates a potential well. The energy of the soliton 
in the presence of the potential becomes 



E = AX^ + ^^+8+^^ (19) 
cosh X cosh X 



When the soliton is far from the center of the potential {X — > oo) (19) reduces to E = + 8. 
It is the energy of a particle with a mass of 8. Some of the features of the soliton behaviour can 
be found from the (19). For example, suppose that a potential barrier of height e is located at 
the origin. A soliton with a low velocity reflects back from the barrier and a high energy soliton 
climbs over the barrier and passes over it. So we have a critical value for the velocity of the soliton 
which separates these two situations. The energy of a soliton in the origin (X=0) comes from 
{\9)E{X = 0) = (4 + 3e) + 8 + 2e. The minimum value of the energy for a soliton in this 
situation is = 8 + 2e. On the other hand, a soliton which comes from the infinity with initial 
velocity Vc has the energy oi E {X = oo) = 4?;^ + 8. It is clear that it can pass though the barrier 
if Vr, > 



2 • 

Equation (18) has an exact solution as follows 

3X^ + 2 _ corili' X (cosh- Xq + f ) 
3Xo' + 2 ~ cosh2^o(cosh2x + f) 

where XO and Xq are initial position and initial velocity respectively. Many of the characters of 

soliton-potcntial system can be extracted from the above solution. In the next sections some results 
are discussed and also compared with the results of the other models. 



4 Comparing of the models 

These three models can be compared numerically. All these models (for a delta-like potential) 
have a parameter in their equation of motion, cr in model 1, Aq in model 2 and e in model 3. The 
parameters control the strength of the external potential. 

It is possible to compare these three parameters in a specific situation by simulation and adjust- 
ing parameters to have same results by different models for that specific situation. It is expected to 
find approximately the relation between the parameters in other situations. A set of simulations for 
the three models have been performed for finding Vc with respect to different values of the potential 
strength with using three models. It is observed that model 3 predicts the value of Vc = when 
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the soliton is far from the center of the potential [Xq —>■ oo). Simulations using models 1 and 2 
show the same behaviour. An effective strength is found by interpolation of simulation results 

on the \J "^2'^ for both models 1 and 2 with respect to parameter of model 3. Figure la shows the 
results of simulations for (1) with sine-Gordon model. Effective strength of model 1 with respect 
to model 3 is 

£1 = (0.0275 ± 0.0022) + (0.786 ± 0.0064)cr (21) 

with standard deviation of 8.5 x 10"^. Figure lb presents the results of the fitting for the model 
2. The result of the fitting is 



62 = (-0.0645 ± 0.00221) + (0.8004 ± 0.0065)Ao 



(22) 



with standard deviation of 6 x 10 ^ 

Figures 1 show that the three models are in agreement with each other if, Vc 



where 



the effective parameters are calculated for models 1 and 2 with respect to parameter of the model 
3. Simulations have been done using Ronge-Kutta method for time derivatives and finite difference 
method for space derivatives. Space grids have been chosen Ax = 0.01, 0.05 and some times 0.025. 
Time cells have been chosen At = in the simulations. Delta function was simulated by the 
function -,/^e~"^^ with several values for a. 
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Figure 1: The critical velocity Vc respect to strength of the potential. Figure la presents the results of 
model 1 and figure lb for model 2. Dotted plots denotes the fitted curve on the function \J^^^ and 
solid lines with data points show the simulation results. 



5 Soliton-barrier system 

A soliton-barrier system is modeled with e > in (18) or (20). Consider a soliton with initial 
velocity of Xq at initial position of Xq — —00. Equation (20) shows that the soliton reaches the 
infinity again with the final velocity X — ±Xq. The soliton goes to —ool+oo) if its initial velocity 
is less (more) than the critical velocity Vc- If the soliton is located at some position like Xq (which 
is not necessary infinity) the critical velocity will not be Neither model 1 nor model 2 has 

analytical prediction for the critical velocity in this situation. However we can investigate this 
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situation numerically using these models. Now let us study the situation with model 3. The soliton 
can pass over the barrier if the soliton energy is greater than the energy of a static soliton at the 
top of the barrier. So a soliton in the initial position Xq with initial velocity of Xq has the critical 
initial velocity if its velocity becomes zero at the top of the barrier X = 0. Consider a soliton 
with initial conditions of Xq and Xq. If we set X = and X — in equation (20) then Vc = Xq. 
Therefore we have from (20) 



/ e cosh^ Xq — I 
2cosh2Xo + % 



(23) 



Figure 2a presents the critical velocity as a function of initial position (Xq) with cr = 0.4,0.6 and 
0.75 in model 1. The equivalent potential strength using model 3 from (21) is ei = 0.3419,0.4991 
and 0.6170 which are shown in the figure 2a with solid lines. Dotted lines in figure 2b show the 
critical velocity as a function of initial position {Xq) using model 2 with Aq — 0.5, 1.0 and 2.0. The 
solid lines in figure 2b show the equivalent situations using model 3 with £2 = 0.2357, 0.7359 and 
1.5363. These figures show that model 3 is in a very good agreement with model 1 and model 2. 




Figure 2: 2a shows the Critical velocity Vc respect to initial position Xq simulated with model 1 and 
figure 2b presents the simulation results for model 2. Solid lines with data points present simulation 
results and dashed lines are plotted using equation (23) with equivalent effective strength calculated 
from (21) and (22). 

If the initial velocity is less than Vc then exists a return point in which the velocity of the soliton 
is zero. This point is derived from (20) 



Xstop = cosh ^ 



'-^ , a = (SX^ + 2) 1 + ) (24) 

2q-4 / ' V J \ 4cosh^Xo/ 



where Xq and Xq are initial position and initial velocity respectively. If the above equation is 
rearranged as 

T = ^+^of- + (25) 

cosh Xstop cosh Xq V ^ 2 cosh Xq ) 

one can find a linear relation between — r^ir; and Xk . Figure 3a shows — as a function 

cosh^ Xstop ^ cosh'^ Xstop 

of Xq with constant value for Xq and some values of e using model 1. All the simulations result the 
same value at = which is equal to ^^^^2 ■ Equation (24) also shows another linear relation 
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between t — ^"^^ cosh^ Xo ' figure 3b demonstrates the numerical simulations with model 1 
for this situation. Model 1 is in agreement with linear relation between — — and — A „ as 

^ cosh'' Xstop cosh'' Xo 

well as linear relation between — t-t^ and Xn, which are conclude from the analytic model 3. 

Model 2 also show the same linear relations. 
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Figure 3: Figure 3a shows linear relation between — — and Xq^ . Linear relation between 

COS II -^stop 

— — and — ,2 ^ has been shown in figure3b. Simulations have been done with model 1. 

cosh Xstop cosh Ao 

The trajectory of a soliton during the interaction by the potential, X(t) follows from (20) as 

-1 



t : 



rxit) / (3x2 + 2) (eosh^Xo + ^) ,„,i^2^ ^1 
Jx{t=o)\\ 3cosh2Xo cosh"X + f 3 I 



The above integral has been evaluated numerically by using Rubmerg's method and X(t) was 
plotted versus t. This result was compared with direct simulation using model 1. Figure 4 shows 
the result for a system with e = 0.4, Xq = —5 and Xq = —0.5. There is a little difference between 
the predicted final velocities from different models after interaction. The difference is reduced when 
the height of the potential(e) reduces. The difference is due to the approximation which is used 
for deriving (15) from (14). Same results have been found when the soliton reflects back after the 
interaction. Another interesting experiment is finding the time that a soliton needs to reach a fixed 
point when it has different initial velocities. This situation has been investigated with both model 
1 and model 3. Figure 5 shows the results for some different soliton-potential systems. 

Several different simulations have been set up and the results of the three models compared. All 
the simulations show very good agreements between three models. Collation between model 2 and 
analytic model 3 also shows a good agreement between these two models, as well as what we can 
see between the model 1 and model 3. In order to avoid adding more figures, the results of model 1 
has been reported. So we can conclude that the analytic model 3 can predict the characteristics of a 
soliton-barrier system. It is concluded that the soliton 'lives' like a point particle but the extended 
nature of the soliton induces some effects on the potential and therefore the effective potential 
becomes more complicated than what we see in a point particle-barrier systems. 
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Figure 4: Trajectory of a soliton during the interaction with a delta-Uke potential. SoUd hne presents 
the result of model 1. Dashed lines show the same situation calculated with analytic model 3. 




Figure 5: The time that a soliton at initial position Xq = —5 needs to reach the final position X=0 
as a function of soliton initial velocity. Solid lines show the simulation results using model 1 and the 
dashed lines present same situation calculated with model 3. 
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6 Soliton-well system 



The soliton-well system is a very interesting problem. The behaviour of a soliton during the 
interaction with a potential well is very different from a point particle in the same situation. It is 
found that some differences can be explained by the characters of the effective potential. 

Changing e to — e in (20) changes potential barrier to potential well. The solution for the system 

is 

3X^ + 2 _ cosh^ X (cosh^ Xq - f ) 
iX^ + 2 ^ cosh^ Xq (cosh^ X - f ) 

Let's examine the validity of (26) by simulating model 1. Models 1 and 2 have similar behaviour. 
Here the simulations are performed using only model 1. 

Consider a potential well with the depth of e. A soliton at the initial position Xq moves toward 
the well with the initial velocity of Xq. It interacts with the potential and reaches a maximum 
distance from the center potential Xmax- The velocity of the soliton at Xmax is zero. Xmax can 
be found from (24) but with e < 0. Figure 6 shows the results of simulations with e = —0.3, —0.2 
and -0.1. The dashed lines show the results of linear fitting on the simulation data with model 1 
and the solid lines presents the results of model 3 with effective potentials from (21). 



■ (k=-0.3) .44 -6.48X (modell) .455-6.87x (modelS) 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 



Figure 6: The dashed lines with data points show results of simulations using model 1 and the solid 
lines present same situation calculated with model 3. The effective potentials has been calculated with 
(21). 

Simulations are in agreement with linear relation between — \ ^ and - . Also there is another 

^ cosh^ Ao e 

linear relation between — ,\ „ and — r^y^ . 

cosn Ao cosn- Xj^ax 

The time required for a soliton with an initial velocity Vq, from initial position Xq to reach the 
origin also has been simulated using model 1 and has been calculated using model 3. Figures 7 
show the results for a soliton with an initial position ATo = —5. The effective potential from (21) 
has been used in model 3. 

These results show that the equation (21) is valid for the potential well too. Also we see that 
the model 3 covers the bahaviour of soliton-well system with an acceptable precision. 
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Figure 7: The time needed for a soliton to travel from initial position Xq = — 5 to origin as a function 
of soliton initial velocity. FigureTa shows the simulation results using model 1 and same situation 
calculated with model 3. The strength of the potential is e = —0.4. Figure 7b presents the results for 
e = -0.3. 

A soliton can pass through the potential well if it has suitable initial velocity. Figure 8 presents 
the simulation results using model 1 and calculation using model 3. 

A soliton with a low initial velocity might get trapped by the potential and oscillate in the well. 
The Period of the oscillation can be calculated with model 3 from (24) with a negative potential 
strength (e < 0). 

Figure 9 demonstrates an oscillating situation. This figure presents the trajectory of a soliton 
in a potential well with a = —0.4 for the model 1 {eeffecUve = —0.29 for model 3). The soliton is 
located at the initial position Xq = —3 and starts moving with an initial velocity Xq = 0.01. The 
period of the oscillation simulated by model 1 is , T « 398 while the period calculated by (27) is 
about 372. 

Due to using adiabatic approximation in model 3, the results show noticeable diff'erences among 
the models when the velocity has rapid or big changes. In a situation where a soliton moves from 
an initial position very far from the potential with a very low velocity the models fail to match by 
using the fitting equations (21) and (22). This means that a better approximation is needed, but 
the analytic model is acceptable. 

An attractive situation in the soliton-well interaction is the fine structure of the islands of 
trapping. In model 1, the final situation of a soliton with an initial velocity lower than the critical 
velocity is very sensitive to its initial conditions and the strength of the well. In most of the cases 
when the incoming velocity of the soliton is smaller than critical velocity the soliton cannot escape 
from the potential. Particularly, after the first interaction soliton will stop and then it will return to 
interact with the potential again. For most of the initial conditions, the soliton will lose its energy 
again during the following interactions and finally becomes trapped by the potential. Equation 
(18) of model 3 clearly shows that the effective force is a function of the soliton velocity, so it may 
be dissipative. However, for some specific initial velocities, the soliton may escape to (±oo) after 
some interactions. Finding an analytical description for the windows of trapping in the model 3 
is hard and we couldn't find a clear analytical formulation for this phenomenon using model 3. 
But this situation can be studied with solving equation (18) numerically. For investigating of this 
phenomenon in the model 3, a plot of "initial velocity respect to final velocity" of the soliton is 
needed. Figure 10 shows the outgoing velocity of the soliton as a function of its incoming velocity. 
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Figure 8: Trajectory of a soliton during the interaction with a potential well. Solid line shows the 
simulation results using model 1 and dashed lines presents the results of the model 3. 
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Figure 9: Oscillation of a soliton in a potential well. Solid line presents the results simulation using 
model land dashed line shows the trajectory of the soliton calculated with model 3. 
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Figure 10: Final soliton velocity as a function of its initial velocity for e = —0.5. Zero final velocity 
means that the soliton is captured by the potential well. 

The initial position of Xq = —2 was used in simulations. The outgoing velocity has been calculated 
when the soliton reaches X = ±10. Zero final velocity means that the soliton is captured by 
the potential. The differential equation (18) has been integrated numerically by using "quality- 
controlled" step size in Runge-Kutta method with maximum error less than 0.001. Simulations at 
this precision show reflection and also transmission from the potential well at some initial velocities. 
We have to check the validity of the results by examining the windows of escaping of figure 10. For 
this purpose we need the trajectory of the soliton with an initial velocity in the region of escaping, 
but with a higher precision. Some simulations with better precisions have been performed using 
Maple which contains some advanced algorithm with higher precision. The results show that for 
most of the escaping windows the soliton trapped in the well if we simulate the situation with 
better approximation. This means that this situation is very sensitive to the precision of numerical 
calculations. However there are some values of initial velocities in which the soliton can transmit 
over the well or may reflect back. This phenomenon needs a deeper investigation. 

7 Conclusion and Remarks 

In this article, an analytical modc;l for the soliton-potential interaction has been presented. It is 
shown that the model has a very close relation with other models in the way that it is possible to 
flt this model over the other models. The model gives the critical velocity in the soliton-potential 
interaction as a function of initial conditions of the soliton and the characters of the potential. The 
model predicts specific relations between some functions of initial conditions and other functions of 
final state of the soliton during the interaction. Also the model presents a good approximation for 
the trajectory of the soliton during the interaction. The oscillation period of the soliton in the well 
can also be calculated by the analytical model. Simulations using other models are in agreement 
with the present analytic model. But this model does not predict the narrow windows of soliton 
refiection from the potential well. 
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The model can be used for prediction the results of other potentials beside the sine-Gordon 
model. 
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